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An analytic solution for an uniaxial spherical resonator is presented using the method of Debye 
potentials. This serves as a starting point for the calculation of whispering gallery modes (WGM) 
in such a resonator. Suitable approximations for the radial functions are discussed in order to best 
characterize WGMs. The characteristic equation and its asymptotic expansion for the anisotropic 
case is also discussed, and an analytic formula with a precision of the order 0[i/~^] is also given. 
Our careful treatment of both boundary conditions and asymptotic expansions makes the present 
work a particularly suitable platform for a quantum theory of whispering gallery resonators. 
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I. INTRODUCTION 

London's St. Paul's Cathedral is famous for its rich 
history and architecture; one of the most unique aspect 
of this building is the whispering gallery that runs along 
the interior wall of its dome . When sounds are uttered 
in low voice against the wall, sound waves generated cir- 
culate around the wall many times before fading away. 
As these waves propagate, they bring with them sounds 
that are audible on the opposite side of the dome. On 
the contrary, if the same sounds are uttered at higher 
volume, the frequencies of these sounds waves will not 
match and a lot of noise is created, making the message 
difficult to be heard at any part of the wall. 

The physical explanation of this effect was firstly given 
more than a century ago in terms of reflection of acous- 
tic rays from a surface near the dome apex. It was ini- 
tially assumed that the rays that propagate along dif- 
ferent large arcs of the dome in a form of a hemisphere 
should concentrate only at the point diametrically oppo- 
site to the source of the sound. Afterwards lord Rayleigh, 
in his Theory of sound [2 , provided a different expla- 
nation of the effect that he named Whispering Gallery 
Waves: sound clutches to the wall surface and creeps 
along it without diverging as fast as during the free space 
propagation: these sound waves then propagate within a 
narrow layer adjacent to the wall surface. It was then 
discovered, at the beginning of the last century, that 
optical whispering gallery waves can exist even in di- 
electric spheres [3J H] . An optical resonator that shows 
this particular wave structure was then called Whispering 
Gallery Resonator (WGR). In recent times whispering 
gallery waves have found new fame with the development 
of nano-optics, in particular with the ability to manufac- 
ture spherical and toroidal WGR with very high quality 
factors that ranges from lO'^ to 10^° [SHE]- This moti- 
vated a large theoretical and experimental work around 
these devices (see for example Ref. IMTT1 and refer- 
ences therein). The ability to store light in microscopic 
spatial volumes for long periods of time (due to the high 
Q-factor) resulted in a significant enhancement of non- 



linear interactions of various kinds like four- wave mixing 
[12 [12] , Raman [14j , parametric and Brillouin scattering 
[TTl [T5] , microwave up-conversion [1^ , second and third 
order harmonic generation [TTj - fTO] . Besides the field of 
nonlinear optics, WGRs were recently used even for cav- 
ity QED experiments (201 HI] ■ For an exhaustive review 
on the applications of WGRs see Ref. [32] . 

Since many of the applications of these resonators in- 
volve nonlinear optics, WGR are commonly fabricated 
using nonlinear materials or anisotropic crystals [llj . 
Despite the wide scientific production in the theory of 
anisotropic spherical resonators, that ranges from gen- 
eralization of scattering methods [531 potential 
method dyadic Green function approach and 
Fourier-based analysis [29], and even though an exten- 
sive study of isotropic WGRs was done in the past [33] . 
detailed studies on anisotropic WGRs are still very few. 
To the knowledge of the authors anisotropic WGRs are 
mainly reported in literature as studied with FDTD mod- 
els [301 EI]) cavity loading [32] and direct solution of 
Maxwell's equation with a surface nonlinear polarization 
as a forcing term [TU] . 

In this work, we intend to develop a suitable analytic 
theory for WGRs, starting from a review of the solu- 
tions of Maxwell's equations in an uniaxial spherical res- 
onator, then presenting and discussing its mode structure 
in the limit of small anisotropy, and finally obtaining the 
spectrum of whispering gallery modes sustained by the 
resonator and their structure, discussing how anisotropy 
influence those modes. A detailed discussion on the ap- 
plication of boundary conditions to this resonator is also 
presented, pointing out how to apply correctly these con- 
ditions and discussing some of their basic features that, 
to the knowledge of the authors, is not present in earlier 
works. It is opinion of the authors that this discussion 
is important in order to better understand the physics 
behind this problem. This is the first main result of this 
work. Finally, we introduce a more accurate approxima- 
tion for the field outside the resonator when the index of 
the Hankel function tends to infinity, as we noticed that 
the commonly used power expansion (as, for example. 
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the one presented in Ref. (33]) does not match the exact 
function completely, i.e. it has an additional phase factor 
respect to the real function. Such phase factor becomes 
relevant when field amplitude, as opposed to field inten- 
sity, turn to be fundamental. This happens, for example, 
when one wants to quantize the electromagnetic field in- 
side the resonator, as required for a proper treatment of 
spontaneous emissions processes. Thus, the present work 
may serve as basis for a quantum theory of WGRs. This 
is the second main result of this work. 

This paper is organized as follows: in section II the 
Debye method of potentials for solving Maxwell's equa- 
tions is briefly presented and then used in section III to 
develop the theory of an anisotropic spherical resonator 
for a dielectric uniaxial sphere. In section IV, Whispering 
Gallery Modes (WGMs) are obtained as limiting case of 
the normal modes of the dielectric sphere with high quan- 
tum numbers and their spectrum is discussed. 



II. DEBYE METHOD OF POTENTIALS 
A. Isotropic Solution 

Before considering the problem of an uniaxial spher- 
ical resonator, it is pedagogical to briefly review the 
method of Debye potentials [S^, largely used to solve 
Maxwell's equations in integrable systems. Let us con- 
sider a monochromatic field with an harmonic time de- 
pendence (i.e. , E{x,t) — i?(a?)e~*'^*) in an isotropic 
sourceless domain ft; Maxwell's equations inside fl can 
be written in the following symmetric form: 



V X £■ = -ikH, 

V X # = ikE, 



(la) 
(lb) 



where k = uj^/e/c is the wavevector in vacuum and and 
e is the dielectric constant inside il. In order to fully de- 
termine the fields, it is necessary to specify their values 
on the domain boundary dfl. Electromagnetic bound- 
aries are usually of two types: perfectly conducting walls 
(the field is zero on dfl), or open systems, where the 
field components inside f2 and the ones outside must 
match on dfl. The boundary, together with symmetry 
considerations, gives a hint on which is the more suitable 
coordinate system to be used to solve the problem (e.g. 
spherical coordinates for spheres, cilindrical coordinates 
for wires etc.) [53] . 

Let us specify our problem by considering an open sys- 
tem constituted by a sphere of dielectric constant e and 
radius R surrounded by an isotropic medium (i.e. air) 
[35] . The set of Maxwell's equations ([T]) in the spherical 
reference frame can be written in the following compact 
form: 



d d 

— — (im£'m) — — - — (LnEn) — ^ik€inmLnLmE[l,{2a) 
OQn oCm 

d d 

-TZl^iLmHm) — TT^ — (LnHn) = ikeinmLnLrriEl, (2b) 

where {l,n,m} € {1,2,3}, Cm are the spherical coordi- 
nates (Ci = ip, (2 = and = r), and are the metric 
coefficients of the spherical reference frame {Li = rsin6', 
L2 — r, L3 — 1). The Levi-Civita symbol on the 
right-hand side of Eqs. ^ is equal to 1 if {l,n,m} is 
equal to {1, 2, 3} or any of its even permutation, is equal 
to —1 for any odd permutation of {1,2,3} and equal to 
zero elsewhere. 

In this reference frame, the fields can be decomposed in 
the so-called transverse electric (TE) and transverse mag- 
netic (TM) waves: TE waves are characterized by having 
Er = 0, i.e. the electric field is transveral with respect to 
the radial direction r. TM waves are instead character- 
ized by having the magnetic field transveral with respect 
to the radial direction (i.e. Hr — 0) [35. For the sake of 
simplicity, let us fix our attention on TM waves; the cal- 
culations for TE waves can be straightforward obtained 
by analogy. From Eqs. (2a) for Z = 3, by substituting 
Hj^ = 
that 



it is possible to introduce the W function such 



r sin 9E^ 



rEo 



dW 

dip ' 
dW 



(3a) 
(3b) 



By substituting relations ^ into Eqs. ( 2b I for Z = 1,2 
and writing W = dU/dr, where U is the TM Debye 
potential, from (2b) we obtain 



IdU 
''r~d0' 



Hfj = ~ik 



1 dU 
r sin 9 dip ' 



and according to Eq (2a I, the r-component (/ 
the electric field is given by 



Er = — 



d f I d 



dp V sin 9 dp 



(4a) 
(4b) 
= 3) of 

U. (5) 



Note that the differential operator that acts on the poten- 
tial U in this equation is the angular momentum operator 

L — —(fx V), that is the same operator that originates 
the centrifugal potential in the Hydrogen atom [3^ . 

Therefore, all the components of the electric field are 
expressed in terms of the U potential solely. In order to 
explicit them, it is necessary to find the equation which 
the U potential satisfy. To do this, we can use one of the 
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last two equations left available from Eq. (2a I, i.e. the 
ones with I = 1,2. By using one of them it is possible 
to obtain the following wave equation that U must satisfy 



decreases to zero as the inverse square of the distance 
for a; 00. Putting everything together, the TM Debye 
potential reads as follows: 



Qj.2 



(6) 



where V^^ is the angular part of the Laplace operator in 
spherical coordinates, i.e. the angular momentum oper- 
ator L. 

The solution can be easily found with the method of sep- 
arations of variables; writing the potential as U{r, 0, f) — 
"^{r)Q{9)^{ip) and substituting this into Eq. (|6|, we ob- 
tain the following equations for the functions ^'(r), Q{9) 
and $((^) 34]: 



'd9 



deJ 
d^ 

dip^ 



7> 

_C2 
. 9 

sm 



= 0, 

e = 0, 



C2$ = 0, 



(7a) 
(7b) 
(7c) 



where Ci, C2 and C3 are the separation constants ap- 
pearing in the equations by separating the variables, 
whose values must be ci — nin A- 1) and C2 = in 
order to the solution to these equations to be unique, i.e. 
physically meaningful; n and ni are integers, including 
zero. With these values the solution is straightforward. 



Equations ( 7b I and ( 7c I give rise to the so-called spher- 
ical harmonics Ynra{0,ip) = ^P^'(cos 6')e""'^, i.e. the 
eigensolutions of the angular momentum operator [53] . 
where P™(cos6') are the associated Legendre functions 
of the first kind that are solution of Eq. (7b), 
the complex exponential is a solution of Eq. 



while 



(7c|. ,yr 



is a normalization constant that guarantees that the 
integral over the solid angle is unitary. The radial 



equation (7a I can be transformed into a Bessel equa- 
tion by the substitution ^'(r) = VkrZ{kr) that brings to: 



f/^/^^*(r,0,^) = Q„t/e,tV^Z,(fcr)y„„(0,^), (9) 

where n, m are the angular quantum numbers that ad- 
dress the single mode of the resonator, Z^{kr) is the 
radial Bessel-type function that is equal to the Bessel 
function J^(kr) inside the dielectric sphere, and is equal 
to the Hankel function of the first kind Hjy^\kor) out- 
side the dielectric sphere. The constants C'int/ext ai'c to 
be determined by applying suitable boundary conditions. 
Note that the argument of the Bessel function inside the 
sphere contains the sphere dielectric constant e via the 
wavevector k = uj^/e/c while the argument of the Han- 
kel function outside the sphere contains only the vacuum 
wavevector ko = ur/c because Sair = 1- 
The components of the electric and magnetic fields for 
TM waves can then be written as a function of U as fol- 
lows: 







(10a) 


Ee 


1 d'^U 
r drdO ' 


(10b) 




1 d^U 
r sin 9 drdip ' 


(10c) 


Hr 


= 0, 


(lOd) 


He 


, IdU 

r oip 


(lOe) 


Hip 


, IdU 
r oO 


(lOf) 



Note that in obtaining the expression of Ej. we have com- 
bined Eqs.([5]) and (|6|. 

If we proceed in a similar manner for TE waves, we ob- 
tain: 



d^Z _^ldZ ^ / )z -0 (8) 

dx"^ X dx \ x"^ ' 

where x = kr and = n + 1/2; the solutions to this 
equation are the four Bessel functions Ju{x), N^{x), 
H'i}\x) = J^{x)+iN^{x) eindHi^\x) = J^{x)-iN^{x). 
Physically, the solution inside the sphere must be finite 
at the origin, and the only plausible solution is Ju{x) 
because N^{x) has a divergence at the origin. Outside 
the sphere, instead, the solution should have the form 
of a runaway wave with the Sommerfeld condition at 
the infinity (i.e. , the solution must drop at infinity as 
the inverse square of the distance). For this reason the 
correct solution in this domain is the Hankel function 
of the first kind Hi^\x) because its asymptotic form 





= 0, 


(11a) 


Eq 


, IdV 

r oip 


(lib) 


E^p 


, IdV 
r oO 


(11c) 


Hj- 




(lid) 


He 


1 d^V 
r drdd ' 


(lie) 


Hip 


1 d^V 
r sin 6 drdip ' 


(llf) 



where V is the TE field potential obtained by Eqs. ^ 
by substituting the TE ansatz Er = 0. 
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B. Boundary Conditions 

Prior to investigate the structure of the modes for 
the anisotropic resonator, it is important to discuss the 
boundary conditions that have to be apphed to this prob- 
lem. At the resonator surface r — R, the wave vector k 
inside the dielectric sphere has to match the wave vector 
ko = uj/c outside the sphere and the constants Cint and 
Cext should be chosen properly. 

There is not a unique way to fulfill boundary conditions: 
in fact, one could apply "pure" or "mixed" conditions: 
the former consists in applying the boundary conditions 
to all the components of only electric or magnetic field, 
while the latter applies the boundary to certain compo- 
nents of one field and certain other components of the 
other field. Obviously, these two different paths bring 
to the same physical solutions [Mj. Among these pos- 
sibilities, in this work we chose to apply "pure" bound- 
ary condition, i.e. we impose that the tangential electric 
(magnetic) field components for TM (TE) waves has to 
be continuous at the resonator surface r — R, while the 
radial component of the displacement vector D = eE is 
continuous across the resonator surface. For TE fields, 
the radial condition is automatically fulfilled, since the 
resonator is non-magnetic (i.e. , /i = 1). 

The condition for the radial component of the dis- 
placement vector {sintE"^*" = SextE^^*) across the 
resonator surface gives the ratio between the inner and 
outer coefhcients 



int 



h\}\kaR) 



(12) 



while the continuity of the tangential component 
= Eg^^ of the field gives rise to the so called 
characteristic equation, that allows to determine the 
allowed values for the wave vector k (i.e. , to find the 
spectrum of the allowed modes) inside the resonator, 
and it turns out to be 



UkR) hi'\koR) ' 



for TM waves and 



[jAkR)]' _ 1 [hi'\koR)]' 



J.ikR) 



re hi^\koR) 



(13) 



(14) 



for TE waves. In these equations jiy{x) = \/xJv{x) and 
h^\x) = ^/xHi^\x) are the Riccati-Bessel functions, 
and the prime indicates the total derivative with respect 
to the argument on which the functions depend, i.e. over 
kR or koR. 

Although formally corrected, as they are these bound- 
ary conditions do not provide a unique solution to the 
determination of the mode patterns in the resonator. 



In order to better understand this not-uniqueness of 
the solution, let us consider the general structure of 
eqs. (13) and (14). Let T{x) be a piecewise function 
defined across an interface, placed at a; = 1, between 
two region of space, such that T{x) — C'^^'^f{x) for 
X < 1 and T{x) = C^'^^ g{x) for x > 1, with f{x) and 
g{x) two arbitrary real valued and regular functions. 
The constants C'^*''^-' are to be determined by the 
boundary conditions and they must be chosen in such a 
way that the following characteristic equation is satisfied: 



f{x) g'{x) 



(15) 



where the apex stands for the derivative of the two func- 
tions with respect to their arguments. Since this is a 



generalization of the characteristic equations ( 13 1 and 



(14), this equation must hold at the interface between 



the two region of space considered, i.e. its validity is lim- 
ited to X = 1. From eq. ( 15 ) it is clear that if we admit 



that the derivatives f'{x) and g'{x) of the functions are 
equal at the separation interface a; = 1, then the func- 
tions themselves will be discontinuous with a jump that 
has the value of 1/a . On the other hand if we now ad- 
mit that the the functions f{x) and g{x) are equal at the 
separation interface x = 1, then their derivatives must 
be discontinuous, and the magnitude of the discontinuity 
is precisely a. 

The first situation corresponds to require that the 
derivative of the function T(x) is continuous at the 
separation interface (i.e. r'(l)+ = T'(l)^, where the 
plus or minus superscript stands for the expression of 
r(x) for X > 1 and x < 1 respectively). This implies 
that C(^VC^*^ = /'(1)/5'(1) and the function T{x) can 
be written as: 



r/(^) 



r(x) = <^ 



m 



X < 1, 



g{x) X > 1. 



It is then clear that taking T'(x) to be continuous at 
the interface results in a discontinuity in the behavior of 
T{x) while passing through x = 1, whose magnitude is 
/'(l)/g'(l), as is depicted in Fig[l] 

Conversely, the second condition on the functions /(x) 
and g{x) implies that the function T{x) to be continuous 
at the separation interface (i.e. , T{1)^ = T(l)^), we 
have Cf'^VC^'-' = /(1)/5(1) and the hmction T{x) has 
the following form: 



0.5 1 1.5 2 

X [a.u.] 

FIG. 1. (color online) The figure shows the behavior of the 
function T{x) and its derivative T'{x) when the condition 
of continuous derivative at the separation interface x = 1 is 
considered. As can be noted, in this case the function shows 
a discontinuity while its derivative is (obviously) continuous. 
For this example we have used f{x) — Ai(a;) and g{x) — e~^, 
and the magnitude of the discontinuity in the function T{x) 
at the separation interface is /'(1)/(?'(1) ~ 0.5457. 
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FIG. 2. (color online) The figure shows the behavior of the 
function T{x) and its derivative T'{x) when the condition 
of continuous function at the separation interface 2: = 1 is 
considered. As can be noted, in this case the function itself 
is (obviously) continuous while its derivative shows a jump 
discontinuity. For this example we have used f{x) = Ai(a;) 
and g{x) — and the magnitude of the discontinuity in the 
derivative T'{x) at the separation interface is f{l)/g{l) ~ 
0.5457. 



' fix) a; < 1, 

/(I) 



3(1) 



g{x) x>l. 



In this case, taking T{x) to be continuous at the interface 
results in a discontinuity in its derivative, whose magni- 
tude is f{l)/g{l), as Fig[2] underlines. 

In both cases, however, it is not possible to make 
the functions and the derivatives both continuous at the 
same time. This fact makes only possible to obtain the 
ratio between the two constants Cmt, Cext and not their 
explicit value: in order to do that, another condition 
must be applied to the problem. This condition depends 
on the particular problem we are dealing on; in scatter- 
ing problems, for example, the incoming field is known, 
and determines the field pattern on the resonator sur- 
face. In this case Cext is known and the ambiguity is 
removed. Another situation in which the ambiguity is 
overcome is by embedding the whole system (resonator 
plus surrounding medium) in an ideal perfectly reflec- 
tive sphere of big, but finite radius Rq, in such a way 
that the boundary conditions at the metallic surface will 
completely determine the fields: this second approach is 
very useful if we are dealing with the quantization of the 
field in such a system. 

We want to end this discussion by pointing out that 
the first situation (derivative continuous at the interface) 



corresponds to the boundary condition for the electric 
field across a dielectric surface: the normal component 
with respect to the separation surface is discontinuous 
by a factor equal to the ratio of the two dielectric con- 
stants of the two regions, while the tangential compo- 
nents (i.e. , the derivative of the radial field in our spher- 
ical case) is continuous at the interface. The second situ- 
ation, instead, corresponds to put continuous the normal 
component of the displacement vector across the separa- 
tion surface, resulting in a discontinuity of the tangential 
component of the displacement vector at the interface. 
While the former correspond to the usual way of impos- 
ing boundary conditions in an electromagnetic problem, 
the latter is never used, but still valid. 

In this work, however, we are neither interested on 
scattering problems nor on field quantization, and so in 
the rest of the paper this ambiguity will not be removed. 
This does not create too much problems because we are 
only interested on the mode structure of the resonator. 
We leave this problem of not-uniqueness to future works. 



III. NORMAL MODES OF AN UNIAXIAL 
SPHERICAL RESONATOR 



Let us consider the same dielectric spherical resonator 
of radius R of the previous section, but with an uniaxial 
anisotropy along the z-axis described by the following 
dielectric tensor: 
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d 




£2:0; (xx + yy) 



(16) 



In order to use this dielectric tensor in Eqs. it should 
be converted in spherical coordinates; this operation 
is simply done by converting the cartesian dyadics 
XX, yy and zz into the spherical dyadics ff , 90 and ipip 
using the standard cartesian-to-spherical transformation 
relations [34]. By performing this transformation, the 
dielectric tensor in spherical coordinates reads 



(17) 




{rHe) 



d 



— {sin Or H^) = 

= —ikor sin6[err{rEr) + EreirEg)]. (19c) 

As in the previous section, we solve the problem for 
the quasi-TM component of the field (i.e. Hr = 0); 
the quasi- TE solution is again obtained using similar 
arguments. We follow the solving procedure described 
in Ref. 37 . Equation (18c) defines the W function as 
in ([3]). Let us combine (19b) and the derivative with 



respect to r of ( 18a 



or sm dip 



with e± = [Szz i £a;x)/2 and e±_ — e^x- We have then 
defined =£++£- cos(26'), egg cos(20) and 

£rg — £-sin(26'). The fact that the tensor components 
depend on the polar coordinate 9 makes the problem to 
find the eigenmodes of the spherical resonator much more 
difficult. Moreover, in an anisotropic system it is in gen- 
eral no longer possible to divide the electric and magnetic 
fields in their TM and TE components. In order to over- 
come the latter problem, we will focus our attention on 
the case of small anisotropy, i.e. A = £_/£+ ^ 1 (this 
approximation is very good if we consider, for example, 
a dielectric sphere made of Lithium Niobate (LiNbOa) 
for which we have Sxx = 5-3, £22 = 6.47 and therefore 
A ~ 0.01). In such a way the fields can be decomposed in 
quasi-TE and quasi-TM oscillations, allowing us to solve 
the problem using the method of Debye potentials illus- 
trated above. To this aim, and for the sake of clearness, 
let us rewrite the set of Eqs. ([2| for the anisotropic case 
as follows: 



— (rsmOE^) - Tri^r) = ikosm9{rHg), (18a) 
or oip 

— {Er) - Q-^rEg) = ikoirH^), (18b) 
Q^irEg) - —irsmOE^) - ifcor sin0(ri/,), (18c) 



and 



lirsineH,)-^iH.) = 

= ikosme[ergirEr) - £gg{rEg)], (19a) 



If we substitute the expression of ^{rHg) obtained from 
the first equation into the second one and if we define 
the differential operator Ix = d^/dr^ + fc§£j_ we obtain: 



dEr 
dr 



IxW. 



(20) 



to r of ( 18b ) gives 



Combining now (19a), and the derivative with respect 



d_ 
dr 



( 1 dW\ 

(r sin QH^) = ik^r sin Q [srdEr - ~'^^^~Q(f) ' 



drde 



(Er) 



dr 



,irEg 



Again, by substituting the expression for ^(rH^) 
obtained from the first equation into the second one, and 
by defining the differential operator Ig — d^/dr^ + k^egg 
we obtain 



kle^gr 



drde 



dW 



(21) 



As in the isotropic case, we want to define as a 
function of the quasi-TM potential [/, in order to fully 
determine the components of the fields as a function of 
the quasi-TM potential U solely. In order to do this, 
let us compare Eqs. (20) and (21). By noting that 
the differential operator Ig commutes with the operator 
{k^Ergr + d^/drd9) that appears in (21 1, is it possible 



to define, after some simple algebra, the function W as 
a function of the quasi-TM potential U as follows: 



d0 



-iko£±{rE^), 



(19b) 



W : 



d_ 
dr 



{ieuy 



(22) 
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This allow us to write the components of the TM electric 
and magnetic field in terms of the quasi-TM potential U 
as follows l36l: 



rE, 
rE. 



TM 

V 

TM 

e 

TM 



TTi m _ 



J 92 

sin 9 drdif 

0, 



rH. 



TM 



sin 6* dip 
d 

See 





(23a) 


-[hu), 


(23b) 








(23c) 




(23d) 


(ieU), 


(23e) 



de 



l-r-) 
drJ 



IccU. (23f) 



Again, if we proceed in a similar manner for the quasi- 
TE waves we obtain: 



Ej^ 
rEr 



rEl^ 



HeV), 



IJeV, 

drdO 

1 9^ 
sin 9 drdip 
0, 

^^(/>), 
sm 9 0(p 

d ' 



(24a) 
(24b) 

(24c) 
(24d) 
(24e) 

(24f) 



where V represents the quasi-TE potential. 



Equations ( 23 1 and ( 24 ) represent the quasi-TM and 
quasi-TE components of electric and magnetic field in- 
side an anisotropic spherical resonator in terms of the 
TM and TE quasi-potentials U and V. 
The next step consists in constructing an equation that 
the quasi-potentials U and V satisfy, whose solutions give 
the mode fields of the resonator. To do that, let us con- 
sider a general electric and magnetic field, whose com- 
ponents are written as the superposition of the quasi-TE 
and quasi-TM oscillations, i.e. Ei = Ef'^' + Ef^ and 
Hi = Hf^^ +Hf^. Substituting this ansatz in Eqs. and 



( 19 1, after some algebra we arrive at a set of two coupled 
equations for the quasi-potentials U and V that reads 

m 



LhV = 2ie-koT 



dU 
dip ' 



LeU = 2ie-koT 



dV 

dp ' 



(25a) 



(25b) 



where we have defined: 



Lh = {Vl + r%)lg - 2e^kl 



dip'^ 



Le — ^0 + ^^^x 
1 



d_ . .d_ 

sm9d9 '''' 90 



T = rL - 2 



1 

sin^ 9 dp'^ ' 



cos f — sm I 



dr^ 



d9dr' 
A2 



7 2 2 

k^r 



lx + 



+ 2e_(l-A)fcg 



dip"^ 



E = cos(26') 



,d^ 

Qj.2 



Vi 



3 1 



dr 



3 -2^1) -(2.) I + 

l-r—\ -2 — 
dr J dpP' ' 

and A = £_/e+ is the anisotropy parameter. 

From Eqs. (25) it is evident that the anisotropy gives 



rise to a coupling between the two quasi-potentials U 
and V\ this coupling is absent in the isotropic case in 
which the two potentials are independent one each other. 
These equations, in fact, contain the isotropic solution in 
the limit of A = 0. Making this substitution in equation 



(25a I and using the definition of the operator L//, we 



obtain: 



(Vi+r2f,)[eC/ = 0. 



(27) 



Because we set A = 0, the two differential operators Ix 
and Iq are equal, since See = e+ — e_ cos 26* = e+ = e 
and = e, i.e. no anisotropy is present anymore. 
We now define IqU = ^ as the isotropic potential, 
and we assume that this potential can be written in a 
separable way, i.e. yl(r, 6*, (p) = 5'(r)y„m(6', 1^9), where 
^nm(^, </?) are the eigensolutions of the angular momen- 
tum operator, whose eigenvalues are —n(n -f 1) (i.e. , 
Viy„m = -~n(n + l)Ynm)- Substitution of this ansatz 
into the previous equation and consequent simplification 
of the angular part then brings to the following radial 
equation: 



(28) 



that is precisely the radial equation ( 7a ) for the isotropic 
potential, whose solutions are the Bessel functions given 
in the previous section. The same procedure applied to 
the quasi-potential V brings to its isotropic counterpart. 
Since we stated that the anisotropy is small (i.e. , A <C 1), 
then we can use the method of separations of variable 
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to solve the coupled equations (25). We then write the 
quasi-potentials as follows: 



p 



p 



Upir)Yn„^{e,ip), (29a) 



v{r, e,^) = J2 Vp(r-, e,^) = Yl Mr)Ynm{o, (2%) 

p p 

where the subscript p stands for the three indexes n,m 
and q on which the quasi-potential depends; the polar 
index n determines the number of field nodes along the 
polar coordinate 0, the azimuthal number m character- 
izes the nodes in the ip direction and, finally, the radial 
index q gives the number of field oscillations along the 
radial direction r that is related with the solution of the 
characteristic equation. Substituting (29) into (25), us- 



ing relationships (4) and (5) of Ref . [37] and equating the 
terms with equal angular part y„,„(0, i^j), we obtain the 
following set of differential equations for the radial com- 
ponents Up{r) and Vp{r) of the quasi-potentials |38) : 



-"■l.ri "n 



-,TM „, 
-'2.ri-2"n-2 



TM 



2e_m/co[&i,n-it'n-i - &2,n+iWn+i], (30a) 



\/ TE 
A(a2,n-2^ 



, TE \ _ 

2 + 0.-i,n+2'"n+2) - 



- 2ATOfco[6l,n_lU„_l - h2,n+lUn+l\- (30b) 

For the sake of clarity, the expresison of the coefficients 
^TM/TE reported in Appendix A. 

These equations require, in general, a numerical ap- 
proach to be solved. However, in the limit of small 
anisotropy, i.e. A <C 1, a solution to Eqs. (30) can be 
searched in terms of power series in the factor A. The 
zeroth order solution brings (as shown before) to the so- 
lution of the isotropic spherical resonator in terms of the 
Riccati-Bessel functions [see Eq. (28)]. The first order 
solution, i.e. the anisotropic correction we are searching 
for, is obtained by neglecting the terms that are propor- 
tional to A^ in Eqs. (30): however, the resulting equa- 



tions contain in the right-hand side a term that is not zero 
(like in the zeroth order solution) but depends on the 
quasi-potentials Un±i and Vn±i. This coupling among 
neighbor radial modes is a signature of the anisotropy, 
that on one hand breaks the azimuthal degeneracy [the 
azimuthal quantum number appears in the definition of 



the coefficients of Eqs. (30)] and on the other hand re- 
sults in a coupling between radial modes. Although this 
coupling results in an impossibility of an analytic solu- 
tion, it can be demonstrated [37] that these terms are of 
the order A^ and at the first order they can be neglected. 
With this argument, Eqs. (30) at the leading order A 
read 



n{n + 1) l+Un{r) = 



(31a) 



(31b) 



where 71 and 72 are the TM and TE (respectively) 
anisotropic factor given by 



1 - 



2m^ 



1-A 



n{n + 1) 
1 - 4™^ 



4?i(rt + l)-3 n{n + l) 



Eqs. (31 ) have the same structure of the radial equation 
for the isotropic case [31]. The only difference is the 
presence of the 7^ terms that modify the arguments of 
the Riccati-Bessel functions and the quasi-potentials can 
be written as 



{x)Ynm,{0,ip) (33a) 

(x)r„„(0,^) (33b) 

where z^{x) corresponds to ju{x) inside the sphere and 
to h''v\x) outside the sphere. Note also that inside the 
sphere, where the anisotropy exists, x = "fir, while out- 
side the sphere x = k^r (the surrounding medium is still 
isotropic). Substituting these expressions in Eqs. (23) 
and (24) we obtain all the components of the electric 



and magnetic fields in an uniaxial anisotropic spherical 
resonator 



Er 

rEg 

rE,„ 



n{n + 1) 

drdO 
1 5^ 
sin 9 drdif 



Z^{'yix)Yn^rni0,V) 



Z„{-1lx)Yn,m{0,^p) 



(34a) 
(34b) 



Hr = 



rHa = - 



ikQe± d 
sin 6' dtp 
d 



86 



z^(7ia;)y„^„(6',v3) (34c) 
(34d) 

Zu{lix)Yn^rn{0-,v) (34e) 



(34f) 
(34g) 



Zy{nix)Ya^mifi, f) 



for quasi-TM fields. Similar expressions can be written 
for the quasi-TE fields by replacing 71 with 72 , exchang- 
ing the role of the electric and magnetic field and setting 
ei_ = 1. 

The characteristic equation can be found by apply- 
ing the boundary conditions and it turns out to be 

iii = jt/koy/e^) 



>(7ifci?) 



hi^\koR) 



(35) 
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for the quasi-TM waves, and 

- [3Al2kR)]' ^ 1 [h\}\koR)]' 

for the quasi-TE waves. 

As can be seen from the previous equations, in the small 
anisotropy regime, the only effect of the anisotropy is a 
rescaling of the radial coordinate; this is in accordance 
with the fact that an uniaxial crystal shows two different 
refractive indexes: one in-plane (y/Exx) and the other 
out-of-plane (y^eH). Different refractive indexes corre- 
spond to different optical paths, and this is exactly re- 
flected in the rescaling effect of the anisotropy onto the 
radial part of the modes of the resonator. Note also that 
at this level of analysis, the anisotropy doesn't affect the 
angular structure {0 and (p) of the modes. Another dif- 
ference respect to the isotropic case is the value of the 
coefficients on the right-hand side of the characteristic 
equations: while the coefficient for the quasi-TE wave is 
analogous to its isotropic counterpart (if we substitute 
the isotropic dielectric constant e with the anisotropy- 
averaged dielectric constant e+), the coefficient for the 
quasi-TM wave reveals the presence of the anisotropy, 
since it is a ratio between the in-plane dielectric con- 
stant and the anisotropy-averaged one. This is not so 
surprising because for a dielectric uniaxial crystal, only 
the TM component suffer directly anisotropy, while the 
TE component does not, because the crystal is magneti- 
cal isotropic. 

IV. WHISPERING GALLERY MODES 
A. Radial functions with large indices 

With the term whispering gallery mode (WGM) are 
commonly addressed the set of modes with a large index 
n; strictly speaking, the real WGMs are only those for 
which it results that n = m and the radial wavefunction 
shows no roots inside the resonator. However, modes 
with indices n ^ m and with q > I, but close to unity, 
have properties that are close to those of WGMs: this 
means that there is no great difference between a "pure" 
WGM and other modes with nearest indices. 

To study such modes, the first thing we have to do 
is to find a suitable approximation of Riccati-Bessel 
functions for large index. This approximation is useful 
either from the numerical (where computing Bessel 
functions of large index is highly time-consuming) or 
analytical (where the approximation gives the possibility 
to work with easier functions that suit better onto the 
problem) point of view. The appropriate approximation, 
however, should be searched bearing in mind that the 
argument of the Bessel function for a WGM near the 
sphere surface is of the order of its index, i.e. v/x ~ 1. 
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FIG. 3. (color online) Comparison between j,^{x) (solid black 
line) and its Airy approximation from Eq ( |37a[ | (dashed red 
line) for — 1000.5. The approximation holds very well up 
to X ~ I/, while for x larger than v (say for x > 1020) it starts 
to fail. 

By introducing the following change of variables 



the Bessel function inside the dielectric resonator can be 
very well approximated by the Airy function of the first 
kind Ai as follows [55] : 

J.(x)^^/2(0'^'Ai(C), (37a) 

i[>(x)]^-y2(^)^^^|[Ai(C)]. (37b) 

The accuracy of this approximation is of the order i'^^; 
if 1/ exceeds 1000, this accuracy is very satisfactory for 
many calculations, this can be seen in Fig|3j where 
Bessel functions of high order are compared with their 
Airy approximation and in Fig|4j where is shown that 
as v grows, the accuracy of the approximation became 
satisfactory. 

For the solution outside the resonator (the Hankel 
function of the first kind) various approximations are 
available. Here we use the following [40j : 

i[i^(tan))-?7)-f ] , 

hi'\x) = firj) =^ {l + 

y'ltanr? 

~-(^r^ + ^Ar-)+OW'']} (38) 

i/V8tan77 2Ata.n^ rjJ ' ) ^ 

where cos 77 = zz/x; if ^ is large enough (heuristically v > 
1000) the imaginary term inside the curly brackets can be 
neglected. The choice of this approximation rather than 
the one presented in Ref. [33 reside in the fact that while 



10 




1000 



3000 



5000 



7000 



1020 1030 
X [a.u.] 



105( 



FIG. 4. (color online) The figure shows the accuracy a as 
a function of the Bessel index v\ the accuracy is defined as 
the ration between the difference of the true function jv{v) 
and its Airy approximation Ai((^) and their sum, i.e. g = 
- Ai(C*)]/[j^(i') + Ai(C*)], where C,* is <" evaluated for 
X = V. As can be seen, as u grows, the approximation becomes 
more precise. 



FIG. 5. (color online) Comparison between the real part of 



h'^u\x) (black solid line), Eq (38 1 (dashed red line) and Eq, 



(31) of Ref.[33^ (dashed blue line) for v — 1000.5. Our approx- 
imation works very well in the region in which the argument 
is grater than the index (i.e. x > 1010), while the approx- 
imation presented in Ref. [33] is out of phase respect to the 
Hankel function. 



the former is very good when the argument of the Hankel 
function is greater than the index (that is precisely the 
case of the outer functions), the latter is not suitable in 
this region, either for being out of phase with respect 
to Hankel function (as shown in Figs. [5] and [6]) or to 
not approximate in the correct way the original function 
(Fig. [7|. Moreover, Fig. [t] shows that the field outside 
the resonator has all the characteristics of an evanescent 
wave, i.e. it decays exponentially as the distance from 
the resonator surface grows. 

In order to justify this evanescent behavior outside 
the resonator, one can directly solve Eq.(|8| in the limit 
r > R (but still close to the resonator surface), where 
the terms x = k^r outside the derivation symbol can be 
substituted with fcoii, leading to the following equation: 
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\ \ 
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/ / \ 
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- Eq. (38) 
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FIG. 6. (color online) Same as Fig. |5]but the comparison is 
(39) made for the imaginary part. 



whose solution is: 



Z,ix)=Coe-'\ (40) 

where 

\l(k^) ^ ^ ^ A{koR)^ ^ 2kQR ■ 

This is, how we are expecting, the expression of an ex- 
ponentially decreasing field, that is in perfect agreement 
with the hypothesis that the field outside the resonator 
is evanescent due to total internal reflection. 



It can be moreover noted that the oscillatory behav- 
ior of the field components outside the resonator (as de- 
piscted in Figs. [5] and |6] for the radial component of the 
electric field) is not in contrast with this hypothesis, since 
it only represent the behavior of the Hankel function as 
r — >■ oo, i.e. it behaves like a runaway wave whose inten- 
sity is decreasing as 1/r^. In the case of WGM, however, 
no radiation will run away towards infinity since the ex- 
ternal field is evanescent, i.e. the radiation is trapped 
inside the WGM and rapidly decreases toward zero when 
the field goes outside the resonator. 
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-- Eq. (38) 
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-- Ref. [33] 
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FIG. 7. (color online) Same as Fig. [sjbut the comparison is 
made for the absolute value; note in this case how the approx- 
imation presented in Ref. [33] completely fails to approximate 
the Hankel function. 



71 



1 dAi(C) _ e_L 
Ai(C) dC " ^ 



1/3 



2.T 



(42) 



the equation for the quasi- TE field can be deduced by 
this one upon changing 71 with 72 and putting = 1. 
In order to find an approximate formula for the solutions 
of this equation, let us firstly analyze the limiting case 
in which — )■ 00; in this case the right-hand side of the 
equation goes to infinity and the only possible solution 
is that Ai(C) = 0, whose solutions are the zeros of the 
Airy function C,q. Let us denote with A^^ the first order 
correction to these roots; expanding both left-hand and 
right-hand size of Eq ( 42 ) in power series with a first 
order accuracy to terms A^^ we can obtain the first 
order correction to the roots Cg, whose expression is 



B. Angular functions with large indices 

For large indices n, the WGM field is concentrated in 
a narrow interval of angles 9 near 9^ = 7r/2; this makes 
possible to approximate the associated Legendre func- 
tions (i.e. the 6'-part of the scalar spherical harmonics) 
with large indices , with Hermite polynomials with small 
indices as follows: 



Ynra{9.,^) 



2'"y7rw! 



(41) 



Detailed calculations for obtaining this result are shown 
in Appendix B. 



AC, 



-(-) 



2\i/3 



(43) 



where 



a 




2(a;2 



(44) 



and Xq is obtained by substituting the value of the 
first zero of the Airy function {C,q — —2.33811) into the 
definition of C, and inverting that relation with respect 
to X. 

Taking into account the definition of C, the eigenvalues 
of the wave numbers for the anisotropic resonator can 
be represented in the following explicit form 



C. Roots of characteristic equations 

The approximations exploited in the previous section 
are very useful in finding an analytical solution to the 
characteristic equation for the eigenfrequencies of the 
resonator; however, due to the anisotropy, some changes 
in the definition of the variables used above must be 



done. First of all, the x appearing in Eqs. (37) and (38) 



has to be different for the inner and outer functions, 
due to the fact that the anisotropy is confined only 
inside the resonator; we can then define x — k^R as the 
outer variable and by consequence the inner variable 
results to he y — '^i^fe^x. Then, the definition of C 
must be changed into C — (^jvf-^^iy — y). After that, 
by substituting Eqs. (37) and (38) into Eq (35), the 



characteristic equation for quasi-TM field gives [H] 



(« 



1/3, 



71 



(C, + AC,) 
1 ■ 



(45) 



Note that because the quantity AC, is complex, the 
eigenvalues of the wave number are also complex. The 
real part of the wave number then determines the 
eigenfrequencies of the mode. Complex eigenfrequencies 
are fully compatible with the open cavity. As can be 



seen from Eq. ( 43 ) , this approximation has an accuracy 
of v^^l'^. more accurate asymptotic expressions that 
allow the calculation of the positions of resonances of 
the modes in an isotropic dielectric spherical resonator 
have been largely studied in literature (see for example 
Ref . [421 - 147] and references therein) and they were given 
with various accuracy with respect to the index v\ in 
Ref. [17] analytic calculations are carried out to the order 
i^~i/'^, in Ref. the eigenfrequencies are calculated 
with an accuracy of 0\v~'^l'^\, while in Ref.[15] the 
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authors give an expression up to the order 0[z/~^/^]. 
Here we report the anisotropic correction of the formula 
found in Ref.[42] that gives the eigenfrequencies with a 
precision of the order of 0[i^~^] 



1 \l/3/ £^ s3/2 



3 / 1 \ l/3 



, c - , 

xp(V-i)c, + o[.-i]}, 



x(46) 



where P = l/(7i£+) for quasi-TM modes and 
P = ej^/(72e+) for quasi-TE modes. 



D. Whispering Gallery Modes 

We now have all the elements for writing the explicit 
expressions for the radial, polar and azimuthal compo- 
nents of the quasi-TM and quasi-TE WGMs. Taking the 



approximations (37), (38) and (41 ), the equations for the 



components of the quasi-TM fields defined in Eqs. (34) 
become 



n{n + \)y/rn .. 



(2) 



y2E^/2\i/6dAi(C) 



sin 6*2'' 



71 «o 



(47a) 



(47b) 



(47c) 



Hr = 0, 

rH, 



Ai(C) X 



Ai(C) 



sin6'2'"w! V2 

X H^(v^a)e-T"'e*"'^, 

2 



d_ 
do 



H^,(vma)e 2 



(47e) 



(47f) 



for the field inside the resonator, while for the field out- 
side the resonator the expressions are the following: 



r^2"'-y/7r?i;! 
X H^(V^a)e-T"'e*"'^ 
^2e+m 



d r 



2™u;! 



f 1 \4f(??) 

V V sin rjJ drj 



do 



rE, 



H^(y^a)e-^" Je""'^, 

im^/^iiko / 1 \ dfjrj) 
sin6'2^-i/2w! Viysinjy/ d-q 

X H^(V^a)e-T"'e™'^, 
sm6'2'"u>! 



(48a) 



(48b) 



X H^,(V^a)e-T« e*"'^, 
ikoe±\/m 



2'^'w! 

_9 r 



fin) X 



Hu,(\/mQ;)e 2" 



Be) 



(48d) 



(48e) 



(48f) 



Similar expressions can be found for quasi-TE WGMs 
by interchanging the roles of the electric and magnetic 
field in the previous expressions, changing 71 with 72, 
putting ej^ = 1 and changing i with —i. 

Figures [8] and |9] show the behavior of the fundamental 
quasi-TM radial (no nodes in radial direction, i.e. <? = 1) 



and polar ( i.e. 



m) WGM component r Er and 



the "first excited" radial (q = 2) and polar (m = n + 1) 
mode for the same component of the quasi-TM field; the 
physical parameters have been set to be e^x = 5.30, e^z = 
6.47 (LiNbOs) and A = 1064 nm. Note that the radial 
component has its maximum very close to the sphere 
surface (dashed vertical line in Fig. [s]), and its position 
shifts on the left, i.e. on the inner part of the resonator as 
the radial number q increases. The polar part, instead, 
is localized around ^ = 7r/2 in its fundamental state and, 
as m becomes smaller than n, the maxima of the polar 
component tent to repel each other from 6 = tt/2. 

In figures 10 to 13 the intensity distribution of the total 
(47d) electric field of a quasi-TM (i.e. Etm = Err + EgO + 
E^ip) is shown, where the components Ef^ {k = r, 9, ip) 



are given by Eqs. (47) for the field inside the resonator 
and Eqs. (48) for the field outside the resonator; f,0 



and ip represent the unit vectors of the spherical basis 
{r,9,ip). In order to obtain the intensity distribution of 
such a field, one has to sum the square modulus of each 
component of the electric field; however, in this particular 
case, the contribution of Eg and E^p is very small and 
localized at the resonator surface, and the total field is, 
with a good level of approximation, fully determined by 
its radial component. The intensity distribution for the 
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FIG. 8. (color online) Radial part of the fundamental (black 
line) and first excited (red line) WGM for the anisotropic 
resonator; the vertical dashed line indicates the position of 
the resonator surface. The fundamental mode has indices 
n = m = 1000 and q = 1, while the first excited mode has 
the same n and m indices but q = 2. 
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FIG. 9. (color online) Polar part of the fundamental 
(black line) and first polar-excited (red line) WGM for the 
anisotropic resonator. The fundamental mode has indices 
n — m = 1000 and q — I, while the first polar-excited mode 
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FIG. 10. (color online) Intensity distribution of the electric 
field of the fundamental quasi-TM WGM. The WGM quan- 
tum numbers are n = m = 1000, g = 1. 



CONCLUSIONS 

In this work, we have developed a classical-optics the- 
ory for an uniaxial spherical whispering gallery resonator. 
We have presented and discussed the mode structure in 
the limit of small anisotropy for such resonator, and ob- 
tained its spectrum. Moreover, we have furnished a thor- 
ough discussion on the boundary conditions and asymp- 
totic expressions for the electromagnetic field in WGRs. 
Our results may be easily generalized to achieve a quan- 
tum theory of WGRs. 



has n = 1000, m — n — I and g = 1. 



magnetic field components of a quasi-TM mode can be 
straightforwardly obtained from Eqs. (47 1 and (48) or 



by nothing that the Hg component of the magnetic field 
has the same intensity distribution as the radial electric 
field component Er and the H^p component, because of 
the presence of the derivative with respect to 6, has the 
same intensity distribution as the one depicted in Fig. 

m 

As the reader can see from these figures, the field is 
nonzero even after the resonator surface {x = 1); this 
is not surprising because in this region the total field is 
evanescent due to the fact that it has been total internal 
refiected by the resonator, i.e. the field is confined in the 
resonator WGM. 
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APPENDIX A: COEFFICIENTS OF EQS. (22) 

Here are reported the explicit expressions of the coef- 
ficients that appear on Eqs. (30). In order to express 



them in a compact form, let us introduce the following 
quantities: 
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FIG. 11. (color online) Same as Fig. 10 but for q = 2; this FIG. 12. (color online) Same as Fig. 10 but for n — m = 1; 



mode represents the first radially excited WGM. 



g± = kle±, 



In 



1 



2n + 2' 

,2;' 



Tn = r Ix- n{n + 1), 



this mode represents the first polar excited WGM. 



- n{n + 1) 



1 + A 



1 - 4to^ 



4n(n + 1) - 3 



+ (1- A).g+r2}/;-2mV(l-A), 
a^l = 2/„/„+i(n - TO + l)(n - TO + 2) X 

"2 _(2n + 3)r-^+7i(n + l) 



3 = 2/„/„_i(n + m)(n + TO - 1) X 



(2n- l)r— + (n + l)(n- 3) 



With these parameters defined, the a and 6s coefficients 



of Eq. ( 30 ) become: 



APPENDIX B: APPROXIMATION OF SCALAR 
SPHERICAL HARMONICS FOR LARGE 
INDICES 

The equation for the 0-part of spherical harmonics is 
the following 



TM _ rj. 



1 - 4m^ 
4n(n + 1) - 3' 



+ 2.g_m^, 



2g^{n - TO + l){n - to + 2)/„/„+iT:„, 



„ = 25_(n + m)(?i + m - l)/„/„_iT„, 



/„(n - m + 1) rlx - {n + 2) 



62,n = + to) Wa; + (n - 1) 



1 



(^sine^^ + n{n + l) 



/ = 0, 



whose solutions are the associated Legendre functions 
/(6») = (cos 61). Since WGMs are located near the 
equator of the resonator, the correspondent functions 
f{6) will be peaked near the angle 9q = 7r/2; in order to 
find an approximate expression for the polar part of the 
spherical harmonics, let us introduce the new variable 
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a = TT /2 — 9: substituting into equation above gives 



d? df 

-—- - tana- 
da^ da 



n{n + 1) - 



/ = o. 



over the real axis is one. This equation gives the ap- 
proximated form of the associated Legendre functions for 
WGMs; substituting it into the definition of the scalar 



spherical harmonics gives exactly Eq (41) 
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FIG. 13. (color online) Same as Fig. [To] but for n — m = 1 
and q = 2\ this mode represents the first radially and polar 
excited WGM. 

We note that, since the modes are localized near 
the equator, a <C 1 and this allow us to expand 
in power series the trigonometric functions that ap- 
pear in the previous equation, i.e. tana ~ a and 
l/cos^a ~ \ ~V o? ■ Substituting in the previous 
equation, writing /(a) — G(a)e" 1^ and performing the 
change of variables a = x/[(to^-|-1/4)^/''] — we obtain 



d^G 
dx^ 



+ 1) -m^ -x^j 



>G = 0, 



Introducing the quantity w = n — m and remembering 
that WGMs are characterized by high values of the 
indices, i.e. n,m ^ 1, the first term that appears 
inside the curly brackets can be simplified as 2w + 1. 
With this substitution the last equation is precisely the 
Hermite-Gauss equation, whose solutions have the form 
G{x) ~ Hyj{x)e~^ Function /(a) then becomes 



/(a) = P;"(cos0) ~ iVi7^,(V^a)e-T" , 

where TV is a normalization factor whose expression could 
be found by requiring that the norm of f{a) integrated 
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